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1 Introduction 

I have known Barry McCoy since 1972, when I spent six months at Stony 
Brook. His enthusiasm for theoretical physics is infectious, and is coupled 
with a great tenacity in tackling long and extremely hard calculations. I am 
still impressed by the evaluation he did then, with Johnson and Krinsky, of 
the correlation length of the eight- vertex model. |I[ In this calculation it is 
not enough merely to find the next-largest eigenvalue of the transfer matrix. 
Instead one has to determine a whole band of eigenvalues, then perform 
the sum over this band. In the thermodynamic limit this sum becomes an 
integration, and one evaluates it by a saddle-point method. A tour de force! 

Much of the work in solvable models stems from the Yang-Baxter relation, 
in particular its star-triangle form. Here I shall show how one can quite easily 
use this relation to show that a certain expression / is invariant, by which I 
mean that it is independent of the rapidities. It is therefore a temperature- 
like variable, being a function only of the modulus and crossing parameter. 
For self-dual Zn invariant models / is one. 

We evaluate / for the Ising, self-dual Potts, self-dual Ashkin- Teller, 
Kashiwara-Miwa, chiral Potts and Fateev-Zamolodchikov models. In all cases 
it is either one, or is related to the modulus k (or k'). 



2 Star-triangle relation 

Consider a model on the square lattice where spins take values 0, ...,N — 
1. Two adjacent spins a and b interact with a Boltzmann weight function 
W pq (a, b) if b is horizontally to the right of a, and W pg (a, b) if b is vertically 
above a, as in the figure. As usual, p and q are rapidity variables, associated 
with the broken lines. 
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To be as general as possible, we should also allow a rapidity-independent 
field weight S(a) for each spin a on the lattice^]. Then the star-triangle 
relation is: 

J2S(d)W qr (a 1 d)W pr (c,d)W pq (d,b) = 1l pqr W pq (c,a)W pr (a,b)W qr (c,b) 

d 

_ _ _ (1) 

J2S(d)W qr (d,a)W pr (d,c)W pq (b,d) = n pqr W pq (a,c)W pr (b,a)W qr (b,c) 

d 

where lZ pqr is some factor independent of the spins a, b, c. 

It is helpful to consider a honeycomb lattice with edge weights W qr , W pr , 
W pq , and a triangular lattice with weights W qr , W pr , W pq , both with site 
weights S. Then the two relations ([!]) are the conditions for them to be 
related, working from either down-pointing or up-pointing stars or triangles. 

For most of the solved two-dimensional models the weight functions W pq (a, b) 
and W pq (a,b) are symmetric, i.e. unchanged by interchanging the spin ar- 
guments a and b. Clearly the two relations ([I]) are then equivalent. The only 
asymmetric model satisfying the star-triangle relation known to the author is 
the chiral Potts model [J, but for this model S(a) = 1 and W pg (a, b), W pq (a, b) 
depend on a, b only via their difference a — b, modulo N. It follows that one 
can convert one of the relations (|l|) into the other by negating (modulo N) 
all spins. 

Thus the two relations (0) seem to be equivalent for all known models. 
We shall not in fact use this, but we shall use a property that is certainly true 
of all the models considered here, namely that there exist spin-independent 
quantities P pq , P pq such that 

uw pq (a,b) = n^(M = p Pq , 

b b 

HW m (a,b) = UW pq (b,a) = P pq , (2) 

b b 

for all values of the spin a. let Y pq be the N hy N matrix with entry W pq (a, b) 
in position (a,b), and Y pq the matrix with corresponding entry W pq (a,b). 

*It may be possible to absorb S(a) into W pq (a, b) and/or W pq (a, b), but to the author 
it seems clearer not to do so. 

^There are other solvable planar models, derived from the three-dimensional Zamolod- 
chikov model, which are asymmetric, but these satisfy the "star-star" relation, rather than 
the star-triangle, pi 
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Then the properties (0) state that the row and column products of Y are 
equal to one another and are the same for all rows and columns. Similarly 
for Y. It is probably not essential to assume these properties: it may be that 
they are implied by ([l]), to within a gauge transformation. However, they 
are in fact satisfied for our specific models, and they do greatly simplify the 
following argument. 

In addition to the matrices Y pq , Y pq , we shall also need the diagonal 
matrices S, X pq \ c , X pq \ c ,X pq \ c , X pq \ c with entries S(a)5(a,b), W pq (c,a)5(a,b), 
W pq (a,c)5(a,b), W pq (c,a)5(a,b), W pq (a,c)5(a,b) in position (a, 6), respec- 
tively. Note that X pq \ c , X' p ,, X pq \ c , X pq \ c depend on the spin c. 

We also define 

8 = [S(0)---S(N- 1)] 1/N . (3) 



3 The factor 7Z pqr 

For the chiral Potts model, Baxter, Perk and Au-Yang |3j conjectured the 
form of the factor TZ pqr . This conjecture was verified by Matveev and Smirnov 
[|]]: their argument generalizes at once to any model satisfying ([[]), as we shall 
now show. 

Regard the spin c as fixed and think of each side of ([]]) as the element 
(a, b) of some matrix. Then the relations take the matrix form 



Y qr SX pr \ c Y pq T^pqr X pq \ c Y pr X qr \ c 

' X Y X' 

qr uyv pr\c 1 pq '^pqr ^-pq\c 1 pr qr\c i 



Y SX i c Y Ttpqr X pq \ c Y X \ c , (4) 



Y being the transpose of Y. Since c takes iV values, (Q) consists of 2N 
separate matrix equations. 

One can now obtain 1Z pqr by taking determinants. Because of (0) we 
obtain not 2N results, but only one, namely 

T^-pqr fqrfpq/ fpr (5) 

(to within an undetermined factor of an iVth root of unity). 
Here f pq is defined by 

/£ = P- 1 det (SY pq ) . (6) 
When S(a) = 1 we regain the formulae of 0, [|]. 
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4 The invariant I 

Define also / by 

T m = Ppl det (SY pq ) . (7) 

We obtained ([|) by holding c fixed and regarding each side of the star- 
triangle relations as elements (a, b) of a matrix product. Suppose instead we 
hold a fixed and think of them as elements (c, b). Then the resulting matrix 
relation is 

Yp r SX q r \ a Y pq T^pqr ^-pq\a ^qr -^-pr\a 

Y pr SX qr \ a Y pq = TZpqr X pq \ a Y qr X pr \ a . (8) 

Taking determinants, these give (for all a and to within an TVth root of unity) 

TZpqr — fpq f pr / f qr ■ (9) 

Finally, holding b fixed and regarding the relations as elements (c, a) of a 
matrix, we get 

V Qlr' — 7? Y' V ~y' 

1 pr "J ^ pq\b 1 qr '^-pqr qr \b 1 pq ^ pr\b 

Yp r SX pq fa Y q r 7Z.pq r Xq r \bYpq X p r \b , (^0) 

TZpqr — fqr f pr / fpq ■ (H) 

The relations (|5]), (H), (|TTD are mutually consistent if and only if 

fpq fpq — fpr fpr ~ f qr f qr ■ (12) 

The only way this can happen is for the quantity 

I = N- 1 f pq -f pq (13) 

to be independent of p and q . (We have introduced the factor 1/N for later 
convenience.) 

Thus it follows, for any model satisfying the star-triangle relations and 
having the properties (|2|), that / is an invariant, independent of the rapidities 
p and q. It is independent of the normalization of W pq (a, b), W pq (a, b), i.e. it 
is unchanged by multiplying them by factors independent of the spins a, b. 
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5 Inversion and rotation relations 

All of the solvable models we shall discuss can be normalized so as to have 
the properties 

W pq (a,b)W qp (a,b) = 1 , W pp (a,b) = 5(a,b)/S(a) . (14) 

In particular, W pp (a, 6) = P pp = f pp = 1. Setting r = p in (JH), then inter- 
changing p with q and using (|5|), it follows that 

YpqSY qp = (fpqfqp) S . (15) 

Further, there always exists mapping R : p —>■ Rp such that 

W pq (a,b) = W q , Rp (a,b) , W pq (a,b) = W q>Rp (b,a) . (16) 

This means that replacing p, q by q, Rp is equivalent to rotating the lattice 
through 90°. For all the models except the chiral Potts model, the mapping 
is very simple: Rp = p + A, where A is a fixed "crossing parameter". These 
relations (|T6| ) imply that 

fpq = fq,Rp ■ (17) 

If K pq is the partition function per site in the large- lattice limit, then it 
follows that 

Kpq ^q,Rp fpqfqp/^qp ■ (1^) 

Together with an appropriate analyticity assumption about fcpq 111 8i domain 
containing the inversion points (p, p) and (p, Rp) , these relations can be used 
(at least for the models other than the chiral Potts model) to obtain n pq 

1, i- 



6 Zjy- invariant models: duality 

Three of the models we shall discuss, namely the Ising, self-dual Potts and 
chiral Potts, are invariant. That is, S(a) = 1 and 

W pq {a,b) = W pq (a-b) , W pq {a,b) = W pq {a-b) , (19) 

the spin-difference functions W pq (i), W pq {i) being periodic in i of period N. 
The matrices Y pq , Y pq are therefore cyclic (Toeplitz). 
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In this case there is a simple duality property |7|, §6.2], || §5]. To within 
boundary conditions, the partition function of the square lattice model is 
unchanged by replacing W pq (j), W pq (j) by the Fourier transforms 



N-l 



n=0 
N—l 

W p d q \j) = N^ 2 £>-*W M (n) , (20) 

n=0 

where uj = exp(2m/N) is the primitive Nth root of unity. 

Using the fact that the determinant of a matrix is the product of its 
eigenvalues, which for the Toeplitz matrices Y pq , Y pq are the above sums, it 
follows that 

PM, Pp q being defined by (@) with W,W replaced by the duals W^ d \ W (d) . 

Making appropriate choices of iVth roots, it follows that / is inverted by 
such a duality transformation, and at the self-dual point 

1 = 1 . (22) 



7 Particular Models 
Ising model 

For the Ising model with horizontal and vertical interaction coefficients J, J, 
N = 2 and 

Y pq ~ I e -K e K J , * pq ~ I e X I ' y Z6 > 

where K = J/ksT, K = J/ksT. Hence 

I 2 = sinh 2K sinh 2K . (24) 



This is indeed the basic invariant of the Ising model - being the modulus 1/k 
or k of its elliptic function parametrization eqns. (2.1a) and (2.1b)], [IDl 
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eqns. (3.67) and (3.71)]. For I 2 > 1 the model is ferromagnetically ordered, 
for I 2 < 1 it is disordered, and when I 2 = 1 it is critical. The Boltzmann 
weights e~ 2K , e~ 2K are functions of the rapidity variables p, q, in fact elliptic 
functions of p — q. We shall not pursue this further as the Ising model is the 
N = 2 case of the Kashiwara-Miwa and chiral Potts models discussed below. 



Self-dual Potts model 

The ordinary Potts model [111 is an N (or q) - state model with weights 
S(o) = 1, 

W(a,b) = e~ K + (l-e- K )5 ab , W(a,b) = + (1 - e-«)5 ab . (25) 

It is solvable at its self-dual point, i.e. when 

(e K -l)(e W -l) = N , (26) 

when it is equivalent to the six- vertex model 0, §12.3]. At this point the 
rapidity parametrization is 

e K . . sin(^ + q-p) ^ = sin(2/x - q + p) 
sin(/i — q + p) ' sin(g — p) 

where /x, the crossing parameter, is defined by 

N 1/2 = 2cos/i . (28) 

For N < 4, fi,p,q are real. For N > 4, they are pure imaginary. As iV — > 4, 
they all become small but enter the equations only via their ratios. 
It follows that 

f P q = e K (l — e~ K ) = 2 cos \x sin(/x + q — p)j sin(2ii — q + p) 

J pq = e K (l - e' K ) = 2cos/x sin(2/i-g + p)/sin(ii + g-p) ,(29) 

and hence 



(30) 



in agreement with 
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Ashkin- Teller model 

The Ashkin- Teller model [|12j is usually formulated in terms of a pair of two- 
valued spins at each site, interacting with neigbouring pairs via Ising and 
four-spin interactions |7], §12.9]. This means that it is a four-state model, 
with 5(a) = 1 and interaction matrices 



Y 



( UJ LOi UJ 2 ^3 \ 

LOi UJq UJ 3 UJ 2 

UJ 2 UJ3 ^0 ^1 

V uj 3 uj 2 Ux uj J 



Y 



I UJ UJl UJ 2 UJ 3 \ 

UJi UJ Uj 3 UJ 2 

UJ 2 UJ3 Ul UJi 

V ^3 ^2 5?1 ^0 J 



(31) 



These matrices are not Toeplitz, but they are cyclic in two-by-two blocks. 
The restrictions @) are obviously satisfied. The weights w Q , . . . ,aJ 3 can be 
expressed in terms of interaction coeficients, but we shall just take them to 
be given parameters. 

By performing a duality transformation on one of the two sets of Ising 
spins, Wegner |L3|], 0, §12.9] showed that this model can be converted to a 
staggered square-lattice eight-vertex model, with weights 



(a,b,c,d) = (ujq + uj\, oo 2 — 0J3, uo 2 + u 3 , uj — io x 
on one sub-lattice, and weights 

(a,b,c,d) = (uo + ui,U2 — uJ3,m — Ui,uJ 2 + uJ 3 



(32) 



(33) 



on the other. 

This staggered model has not been solved and we know of no star-triangle 
relation relevant to it. However, it can be solved when the two sets of eight- 
vertex weights are proportional to one another, as it is then the regular 
model: if 

(a,b,c,d)=£(a,b,c,d) , (34) 

then a, 6, c, d can be parametrized as elliptic functions of p — q so as to satisfy 
([!]). However, we do not need this parametrization to note from ([JT]) - (^) 
that 

detF p9 = 16£ 4 P W , dety M = 16 , (35) 

and hence ( to within choices of Nth. roots) f pq = 2£, / = 2/£, and 



(36) 
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Thus I = 1, as for the self-dual Z/v-invariant models. In fact, this solvable 
case of the Ashkin- Teller model is self-dual [13, 14], but it is not critical. 
We know that the associated regular eight-vertex model can be naturally 
parametrized in terms of elliptic functions of some modulus k flTJ, [16[]. Only 



when this k is ±1 , or oo can the model be critical. 

The explanation of this apparent contradiction is that the general Ashkin- 
Teller model has two critical temperatures, on either side of the self-dual 
temperature defined by (|34|), which map into one another by duality [13|, [14] . 
(The duality relation is easily obtained from the above remarks: two Ashkin- 
Teller models are dual to one another if they transform into staggered eight- 
vertex models differing only in interchanging the two sub-lattices. This means 
that the cj-weights of one are proportional to the eigenvalues of the Y, Y 
matrices of the other.) 

Indeed, when the four-spin interaction vanishes the Ashkin Teller-model 
becomes two independent Ising models. The solvability condition ([34]) then 
implies that these two models are dual to one another, and / is the product 
of the individual Is given by (f24|). This product is necessarily one. 



Kashiwara-Miwa model 



In 1986 Kashiwara and Miwa |T^J presented an N— state generalization of the 
Fateev-Zamolodchikov model that breaks the Z^— symmetry, but retains the 
reflection symmetry and the rapidity difference property. This model has also 
been studied by Hasegawa and Yamada [18|, and by Gaudin [19]. Let K, K' 



be the complete elliptic integrals of the first kind of moduli k, k', and let 
q = exp(—TrK'/K) be the corresponding "nome". Then the elliptic theta 
functions H(u), Q(u) of argument u and modulus k, as defined in section 
8.181 of [ECU, are 



H{u) = 2q 1/4 sm(nu/2K) - 2q 2n C os(ttu/ K) +q 4n )(l-q 



n=l 



Q(u) = H(l- 2q 2n - 1 cos(Ttu/K) + q 4n ~ 2 )(l - q 



(37) 



71=1 



We also use the function H 2 (u) = H(u)Q(u). One can write K, K', k, k! 
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in terms of infinite products involving q, in particular 

Let ( be some arbitrary integer and define two functions 

r(n) fi H[KW-2 + v)IN\ f (n) ft Q[K(2j -2 + v)/N] 
Tv{n) ~ } = \ H[K(2j-v)/N] ' tv[n) ~ /i e[K(2j-v)/N] ' 

(39) 

where v is real and the argument n is an integer. These functions are periodic: 
r v {n) = r v (n + N), t v (n) = t v (n + N), 

Then the weights of the Kashiwara-Miwa model are (for integers a, b ) 

S{a) = Q[2K(2a + C)/N}/Q{0) , (40) 

W pq (a, b) = ri_ g+p (a - b) ti- q+p (a + b + () , (41) 

W pq (a, b) = r 9 _ p (a - b) t q - p (a + b + () . 

Note that for this model, unlike all the others we consider, the site weight 
function S(a) is not unity. Curiously, W pq (a, b) and W pq (a, b) are products 
of functions of the spin difference a — b and of the spin sum a + b. Because 
r v (n) = r v (—n) is an even function of n, W pq (a,b) and W pq (a, b) are each 
unchanged by interchanging the spins a, b. The model is therefore reflection- 
symmetric, i.e. non-chiral. 
Define also 

9 = TT \r ■ (42) 

m=l ^ 1 « 

Then, based on numerical calculations, the known N = 2 Ising case and the 
Fateev-Zamolodchikov result (|6^) below, I conjecture that 

Uq 9 }± H 2 [K(2 J+P -q)/N] ' {A6} 

where N' = [N/2] is the integer part of N/2. (The constant g is determined 
by the requirement that f pp — 1.) 

Since / p9 is obtained from f pq by replacing q — p by 1 + p — g, it then 
follows that 

I = g 2 = k' N K N /(k'K) , (44) 
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where k' N , Kn are the k', K corresponding to the elliptic nome q . This 
explicit expression for / is indeed independent of the rapidities p, q. Note 
also that f pq is a single- valued meromorphic function of q—p, despite the iVth 
root implied by the definition (||): this is consistent with the observation [|J 
that f pq is defined by (§) and (|l|) as a single- valued function of the Boltzmann 
weights (with no Nth root), to within factors that are functions of either p 
or q only. 

The factor f pq f qp in fll^D also simplifies. Define the elliptic function 

H 2 {2Ku/k) ~ l-2g m cos2 M + g 2m 

G(u,q) = — — — — = tanw r^— ■ 45 

V H ' H 2 {K + 2Ku/tt) 1 + 2<? m cos 2u + q 2m K ' 

Then 

f P J qp = Ng 2 G(w,q)/G(Nw,q N ) , (46) 
where w = ir(q — p)/(2N) (c.f. equation |15| above and proposition 2 of Ref. 



0) 



Chiral Potts model 

In 1987 Au-Yang, McCoy, Perk and others found solutions of the star-triangle 



relations for the three- and four-state chiral Potts model |2T], |22j . These were 
generalized to an arbitrary number N of states in 1988 0. 

In this model, the rapidity is a point (a p , b p , c p , d p ) on the homogeneous 
curve H 

a N p +k'b^ = fc< , k'a N p +b N p = kc N p , (47) 
where k and k! are two given constants, related by 

k 2 + k' 2 = 1 . (48) 



Related quantities |23[ are x p , y p , u. p , t p , defined by 



x 



d p / d p , y p b p / Cp , fXp dp/ Cp , (49) 



^"pVp ®pbp/ ' c p d p 

They satisfy 

< + < = fc(l + <yf) , fccf = 1-^ , fcyf = l-fc'< • (50) 
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The chiral Potts model is Zjv-invariant, so we can write the Boltzmann 
weight functions as in ([19]). Setting u = e 2m/N , the functions W pq {ct - b), 
W pq (a — b) are 

(51) 

These Boltzmann weights are positive real if x p ,x q , y p ,y q ,ux p all lie on the 
unit circle and are arranged sequentially in the widdershins direction. 



The function f pq has been evaluated |24|, ^5[] . In the present normalization, 
with W pq (0) = W pq (0) = 1, it is 

f N = TT [ ^ ~ ~ g^gXgg ~ uj yp) V (52) 



/xp(x p - uix q )(y p - uj j y q )(x p - w*y q ) 



When iV = 2 it reduces to the Ising model. 

Let R be the automorphism that takes the point (a p , b p , c p , d p ) to (or p , 
bap, cr p , (Ir p ) = (bp, coa p , d p , c p ). Then the rotation relations (0) are satisfied. 
Using (|T3"D, ( p~T| ) and (|52|), we deduce that 



I = f P J Pq /N = l/k'^ N - 1)/N . (53) 

Thus / is indeed rapidity-independent for the chiral Potts model, being sim- 
ply a power of k! . It is inverted by the duality relation k 1 — > 1/k' || §5]. 

The result (^) is not new, being given in equation (2.47) of Ref. |35f| . 
What we have done here is show that its structure is a direct consequence of 
the star-triangle relation, and is shared by other solvable models. 

We can also deduce that 

f , N($ -t%)(x p -x q )(y p -y q ) 



in agreement with (2.48) of Ref. |25[ 



It is not obvious from (p2|) that the rhs is the iVth power of a single- valued 
function (to within factors that depend only on p, or only on q). However, 
the only singularities whose location depends on both p and q occur at the 
zeros and poles of the bracketted expressions. These can only occur at certain 
points V on the (x q , y q ) or (x p , y p ) Riemann surface, namely when x^ = x^ 
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and y q = i/p, or when x q = y p and y^ = x p . We can count these zeros 



and poles by postulating the existence of functions 0^, 0^ of p and q |p6 
such that 0jj has simple zeros when 

x q = oj' L y p and y q = uj r x p , (55) 
and Qij has simple zeros when 

x q = uj % x p and y q = uj 3 y p . (56) 
Then, to within factors that are analytic and non-zero at all such points V 

N-l N-l 

x q - uj r y p = J[ i>m , y q - uj ] x p = JJ Q md 

m=0 m=0 
N-l N-l 

x q - lo 1 x p = J[ i>m , y q - u 3 y p = JJ Q md (57) 

m=0 m=0 
N-l 

tq U) t p j J ] m,i—mQ ^m.i—m 

m=0 

Formally substituting these expressions into fl5"2|), we find that each 0^- or 
Qij function occurs with a power N, or —N on the rhs, so we can take the 
iVth root and obtain 

N-l N-i 

fm = II U^N^+j/Qij , (58) 
i=i j=i 

to within factors that are analytic and non-zero at all points V. So f pq 
certainly has no branch point singularities at the points V: at worst it has 
simple poles or zeros. 

Fateev-Zamolodchikov model 

When k = we can choose the parameters of the chiral Potts model so that 
x p = e^ N , y p = e *^^ N , fi p = 1 , (59) 
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where p is the rapidity variable. Making the corresponding choices for 
x q ,y q ,u. q , (0) becomes 

W (n) = IT M-( P - q + 2j-l)/2N} 

Wpq[U) /i sin[7r(g-p + 2j-l)/2iV] ' lbUj 

_ » sin[7r(g-p + 2j-2)/2JVl 
/i sin[ 7 r(p-g + 2 J )/2iV] " 

We obtain precisely the same result if we set k = (and hence q = 0) in the 
Kashiwara-Miwa model, so the two models then coincide. 



The resulting model is the Fateev-Zamolodchikov model p7 |. Its weights 
( |S0D are real and positive if < q — p < 1. It is a Z N model, self-dual and 
critical. When iV = 3 it is equivalent to the self-dual Potts model (with 
H = 7i /Q and p, q scaled by //). 

From both (^4|) and (|53|) we deduce that 

7=1, (61) 

in agreement with (p2|). 

We find that the right-hand side of the chiral Potts result (|52|) is now 
indeed a perfect iVth power, and that 

_ /2 " sin[7r( g - P + 2j-l)/2iV] 
/w " /i sin[7r(p-g + 2j)/2iV] ' 

where iV' = [N/2]. This agrees with, and is part of the evidence for, the 
conjecture (f43D for the Kashiwara-Miwa model. 



8 Summary 

The method of Matveev and Smirnov can be used to show that the quan- 
tity 7 defined by ([]), (0), ( |TB| ) is independent of the rapidities p and q. For 
particular models this is not necessarily a new observation, but the realisa- 
tion that it is a direct consequence of the star-triangle relation does provide 
a unifying feature for the known solvable edge-interaction models. We hope 
this is a small step forward in the subject to which Barry McCoy has made 
so many outstanding contributions. 
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